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We present numerical evidence for the approximate 5*0(5) symmetry of the Hubbard model on a 10 site cluster. 
Various dynamic correlation functions involving the tt operators, the generators of the 5*0(5) algebra, are 
studied using exact diagonalisation, and are found to possess sharp collective peaks. Our numerical results also 
lend support on the interpretation of the recent resonant neutron scattering peaks in the YBCO superconductors 
in terms of the Goldstone modes of the spontaneously broken SO (5) symmetry. 



PACS numbers: 7.10.Bg, 71.55.Jv 

Although the single band Hubbard model has been exten- 
sively studied in recent years in connection with high-T c su- 
perconductivity, its low-energy content has so far eluded both 
analytical and numerical investigations. While being decep- 
tively simple, this model may have many possible competing 
ground states, and it has proved to be very difficult to organize 
the low-energy degrees of freedom. Recently, a new analytic 
approach based on a symmetry principle has been suggested. 
It was noticed that the Hubbard model enjoys an approximate 
50(5) symmetry, unifying antiferromagnetism (AF) with d- 
wave superconductivity (SC) [[[[]. This symmetry principle 
gives a simple description of the transition from an AF ground 
state to a SC ground state as the chemical potential is varied, 
and it gives a unified treatment of the low-energy collective 
degrees of freedom of the Hubbard model. 

Motivated by its potential importance to the high-T c prob- 
lem, we undertake a numerical finite size study to test this ap- 
proximate symmetry in the single band Hubbard model. The 
central object of our study is the so-called tt operator defined 
as follows, 



's/d 



^ (cos k x ± cos k y ) 



k 



k+Q,T 



-k,T 



Q = i>,7T) (1) 



This operator carries spin 1, charge 2 and total momentum 
(ir, 7r). Charge conjugation and spin lowering operation gives 
five other similar operators. Here the subscript a = s, d refers 
to the internal s-wave or c?-wave symmetry of this composite 
operator. The tt operators were first introduced by Demler and 
Zhang [|[| to explain the resonant neutron scattering peaks in 
the YBCO superconductors and are constructed by fol- 
lowing the analogy with the rj operators considered by Yang 
[Q] . More recently, it was found [ l[j that together with the total 
spin and charge operators, the six tt operators form the gen- 
erators of the 50(5) algebra, and furthermore, they rotate AF 
and SC order parameters into each other. In reference [||] and 
Hp, it was argued that the Wd operators are approximate eigen- 
operators of the Hubbard Hamiltonian. This implies that the 
50(5) symmetry is an approximate symmetry of the Hubbard 
model [|l]], and one can use it to constrain the form of the low- 
energy effective Hamiltonian. 



In this work, we present exact numerical diagonalisation 
studies of the dynamic correlation functions involving the tt 
operators and the AF and SC order parameters. We verify that 
the tt operators are approximate eigen-operators of the Hub- 
bard model, and show that the properties of the various mixed 
correlation functions are consistent with the anticipated pat- 
tern of the 50(5) symmetry breaking. Our exact numerical 
calculation is performed on a a/TO x \/T0 site Hubbard clus- 
ter, using standard exact diagonalisation (ED) methods based 
on the Lanczos algorithm Mixed dynamical correlation 
functions are calculated utilizing a technique previously used 
to obtain anomalous Green's functions [H]. 

As a check to our numerical algorithm and the general 
framework, we first performed our calculations on the cor- 
relation functions involving the i] operator [Q], its mixed cor- 
relation functions with the charge-density-wave (CDW) order 
parameter |Q] and the dynamic auto-correlation function of the 
CDW order parameter for the negative U Hubbard model. The 
relationship between the rj and the tt operators is explained in 
reference [Q]. Since the 77 operator is a exact eigen-operator 
and the negative U Hubbard model is known to have a s-wave 
SC ground state away from half-filling, various exact results 
can be derived for these correlation functions [^]. Our numer- 
ical results are consistent with these exact results. However, 
due to the space limitations, we will not present these calcula- 
tions here. 

Let us first consider the auto-correlation function of the tt 
operators, defined as follows, 
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Here H is the standard Hubbard Hamiltonian, 
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its ground 



state with N electrons and Eq the corresponding ground state 
energy. 3 takes the imaginary part of a correlation func- 
tion. Throughout this paper, we measure the energy of a spec- 
tral function from the ground state energy of the intermediate 
state. Since 7r]^ d changes particle number by two, the natu- 
ral energy scale of the intermediate states is the ground state 
energy of a N + 2 electron system. 

Figures |[a), 0(b), ||a) and |(b) plot the tt+(uj) correla- 
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tion function for U = 4t and U = 84, with electron densi- 
ties (n) = 0.6 and 0.8 respectively. We see that these dy- 
namic correlation functions are dominated by a sharp peak, 
well separated from a higher energy continuum. If ir d was an 
exact eigen-operator of the Hubbard Hamiltonian, its dynami- 
cal auto-correlation function should consists of a single peak. 
Here we see that it is only an approximate eigen-operator in 
the sense that there is also a high-energy continuum in addi- 
tion to the resonance peak. The ratio of the spectral weight 
under the peak to the total spectral weight is found to be 
0.63, 0.68, 0.46 and 0.31 in figs |l](a) to |(b). The separation 
between the peak and the continuum, and the large relative 
spectral weight of the peak demonstrate that the ir d operator 
is an eigen-operator of the Hubbard model to a good degree of 
approximation. The energies of these peaks are 1.064, 0.664, 
0.554 and 0.244, respectively, and are low compared to U. 
This can easily be understood from the fact that the tt opera- 
tors are spin triplet operators, therefore, the bare U interaction 
should not enter its energy scale. 

In reference [|[, the dynamic auto-correlation function of 
the ltd operators was calculated within the 4-matrix approx- 
imation in the particle-particle channel. Within this approx- 
imation, the dynamic auto-correlation function consists of a 
single peak, with a spectral weight proportional to 1 — n, i.e. 
the hole density, and an energy of J(l — n)/2 — 2/i, where 
J w 44 2 /U is the AF exchange constant and [i is the chemi- 
cal potential measured from half-filling. Strictly speaking, the 
4-matrix approximation should be valid only in the low elec- 
tron density limit. Here we see that the qualitative features 
of the 4-matrix calculation are still valid near half-filling. Our 
exact numerical calculations show that the 7rj"(o;) correlation 
function is indeed dominated by a single peak whose spectral 
weight is proportional to the hole density 1 — n. Furthermore, 
the energy of the n resonance peak is also proportional to the 
hole density, with a scale comparable to J. 

The 7r resonance is a composite particle made out of two 
electrons. Near half-filling, it is very difficult for a single 
electron or hole to propagate coherently. In fact it is not clear 
if there is a quasi-particle peak in the one particle Green's 
function at all. In view of this fact, it is rather surprising 
that a two electron Green's function 7rj~(w) has a coherent 
peak while the individual electrons propagate incoherently. To 
demonstrate that the it peak is a genuine collective behavior, 
we plot in figures |l](c) and ||(c) the bubble approximation to 
the tt^(uj) correlation function. This approximation to Trt (u;) 
consists of a particle-particle bubble with fully dressed one 
particle Green's function of the \/T0 x %/l0 site cluster in- 
side the bubble. It fully takes into account the single-particle 
dressing effects, only the vertex correction is neglected. From 
these figures, we see that the bubble approximation yields a 
broad spectral distribution, without any identifiable resonance 
peak. The height of the broad spectral distribution is one or- 
der of magnitude less than the height of the ir peak found in 
the full calculation. This calculation demonstrates unambigu- 
ously that the collective behavior, or in other words, the vertex 
correction, is responsible for the existence of the ir resonance. 
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FIG. 1. Dynamic correlation functions of the V10 x VlO Hub- 
bard model with U = 4t: (a) vtJ(oj) spectrum at (n) = 0.6, (b) 
7rj"(u;) spectrum at (n) = 0.8, (c) bubble approximation to 7rj~(u;) 
as shown in (b), (d) 7r+(ti;) at (n) — 0.8, (e) off-diagonal correlation 
function Sd(u>) at (n) = (0.8/1.0), as defined in the text, (f) spin 
correlation function xq(^) and (g) bubble approximation xq(<^)- 

Figures [l|(d) and ||(d) give results for 7T+ (w) at U = 44 and 
U = 84, respectively, at density n = 0.8. For U — 84, the 
spectrum is broadly distributed and no resonance peak can be 
identified. For U = 44, although there is an isolated spectral 
feature, it is at a much higher energy (« 24) compared to the 
7rj~(w) counterpart, and the height of the peak is almost one 
order of magnitude less. We interpret these results as evidence 
that in contrast to ir a , t:\ is not an approximate eigen-operator 
of the Hubbard model near half-filling. This result is consis- 
tent with the conclusion of the 4-matrix calculation [[|]. The 
crucial difference between the 7rj"(o;) and nf(uj) correlation 
function has recently been used by Zhang [[!]] to argue that 
there is an approximate symmetry between AF and d-wave 
SC, but no symmetry between AF and extended s-wave SC 
near half-filling. 

Originally, the 7r operators were introduced [[|] in connec- 
tion with the resonant neutron scattering peaks recently found 
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in the YBCO superconductors [ph. Below the superconduct- 
ing transition temperature T c , a resonance peak appears in 
the polarized neutron scattering amplitude with momentum 
transfer (tt, tt), and a resonance energy of 25meV, 33meV 
and AlmeV, depending on doping. Demler and Zhang [|[] 
identified this experimental feature with the tt resonance, and 
showed that the particle-particle tt resonance can be mixed 
into the particle-hole channel below T c , and therefore appear 
in the neutron scattering cross-section. A number of other the- 
oretical papers [^] explain the experimental feature in terms of 
a particle-hole threshold behavior near the superconducting 
gap or possible excitonic states inside the gap. More recently, 
the tt resonance has been interpreted more broadly by Zhang 
[|]] as the Goldstone boson associated with the spontaneous 
breaking of the SO (5) symmetry below T c . 




FIG. 2. Same as in Fig. but here for U 
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In this work, we test the ideas of references [g] and [[!]] by 
studying the mixed correlation functions involving the spin- 
density-wave order parameter 



k 



C k+Q,| Ck <l 



Q 



(3) 



and the TTd operator, defined by 
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and the dynamical spin correlation function itself, defined by 
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The correlation function Sd(w) satisfies an important exact 
sum rule 



/oo 
duS d {u) = 2 «- 2 |A d K) 
-OO 



(6) 



where = X)k( cos ^ — cos ^y) c k.T c -k,i is the d-wave su- 
perconducting order parameter. This sum rule follows from 
the fact that the TTd operator is a 50(5) symmetry generator 
which rotates the AF order parameter into the d-wave SC or- 
der parameter Mixed correlation functions like S d (u>), 
which involve a symmetry generator and an order parame- 
ter, are commonly used to prove the Goldstone theorem on 
the Goldstone bosons associated with spontaneous symmetry 
breaking [^|. 

In figures 0(e) and |(e), S d (u) is plotted for U = At and 
U = 8t, for electron density between n = 0.8 and n = 1. We 
see that the spectrum is dominated by a single peak with rela- 
tive weight 0.77 and 0.72 respectively, located at the same en- 
ergy as the tt resonance. We also verified that the sum rule (^|) 
is satisfied by explicitly computing the d-wave order parame- 
ter defined on the right hand side. Therefore, our result shows 
that the 10 site Hubbard cluster at densities between n = 0.8 
and n = 1 possesses considerable d-wave SC fluctuations, 
and furthermore, there is indeed a Goldstone pole associated 
with this spontaneous symmetry breaking, consistent with the 
SO (5) theory /g/. The finite energy of the Goldstone pole 
results from the fact that the 5*0(5) is not only broken sponta- 
neously by the SC order parameter, but also broken explicitly 
by the chemical potential [[j]]. (Sometimes, a Goldstone boson 
associated with a symmetry which is broken both explicitly 
and spontaneously is called a pseudo Goldstone boson. The 
pion in particle physics is such an example.) We also calcu- 
lated Sd{u>) between the densities n — 0.6 and n — 0.8 and 
found it to be nearly zero. Similarly, the mixed correlation 
function S s (oj) involving the tt s and Sq operators vanishes 
for all densities. This result is consistent with other numer- 
ical evidence for d-wave pairing fluctuations in the Hubbard 
model [0. 

Finally we show our calculation for the spin correlation 
function Xq( w ) m figures [j](f) and ||(f), for U — At and 
U = 8t respectively at density n = 1. We see that there is 
a sharp resonance feature at the same energy as the tt res- 
onance. The fact that all three correlation functions 7rj~(o>), 
Sd{oj) and Xq( w ) nave resonance peaks at the same energy 
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clearly demonstrates that the peak in the spin correlation func- 
tion Xq( w ) i s due to a particle-particle intermediate state 
ir,\$>Q~ 2 ). Because of the finiteness of the mixed correla- 
tion function 5<z(w) and the right-hand-side of equation (^), a 
particle-particle excitation at density n = 0.8 makes a finite 
contribution to the spin correlation function at n = 1. This 
feature confirms the argument of reference [||], which explains 
the experimental resonant neutron scattering cross section in 
terms of the collective particle-particle tt resonance. 

To exhibit the collective nature of the resonance peak in 
Xq(tj), and to contrast it with single-particle behavior, we 
plot in figures HVg) and&g) the bubble approximation to 
Xq(^) on the VlO x \/l0 site cluster. In contrast to the full 
spin correlation function, the bubble approximation does not 
posses any low-energy collective poles. These results show 
that it is necessary to include the vertex corrections and the 
mixing into the particle-particle channels to calculate Xq( w )> 
in order to exhibit its resonance peak. These calculations have 
been performed recently and will be published elsewhere Jill ] . 

In conclusion, we have numerically verified one of the most 
fundamental assumption of the 5*0(5) theory of the Hubbard 
model [jl[], which asserts that the 50(5) symmetry genera- 
tors, the TTd operators, are approximate eigen-operators of the 
Hubbard model both near and away from half-filling. Their 
dynamic auto-correlation function shows well separated res- 
onance peaks with large relative spectral weight and low- 
energy. This behavior is intrinsically collective, and can not be 
reproduced by any calculations which neglects vertex correc- 
tions. In contrast to the TTd operators, the tt s operators do not 
show well separated peaks near half-filling. This result shows 
that there is only an approximate 50(5) symmetry between 
the AF and the d-wave SC order parameters, no symmetry 
between the AF and the s-wave SC order parameters. Close 
to half-filling, between the densities n = 0.8 and n = 1, 
the 10 site Hubbard cluster has finite <i-wave superconducting 
fluctuations, which lead to a non-vanishing mixed correlation 
function involving the 50(5) symmetry generator and the AF 
order parameter. Because of the finiteness of this mixed corre- 
lation function, there is a contribution to the dynamic spin cor- 
relation function from the particle-particle tt resonance. These 
observations are consistent with the proposed explanation of 
the resonant neutron scattering peaks in the YBCO supercon- 
ductors [0]. 
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